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ABSTRACT: The unsteady, compressible, thin-layer Navier-Stokes equations are used to obtain 
three-dimensional, asymmetric, vortex-flow solutions around cones and cone-cylinder configura- 
tions. The equations are solved using an implicit, upwind, flux-difference splitting, finite-volume 
scheme. The computational applications cover asymmetric flows around a 5° semi-apex angle cone 
of unit length at various Reynolds number. Next, a cylindrical afterbody of various length is added 
to the conical forebody to study the effect of the length of cylindrical afterbody on the flow asymme- 
try. All the asymmetric flow solutions are obtained by using a short-duration side-slip disturbance. 

1. INTRODUCTION 

The problem of asymmetric vortex-flow around slender bodies has received considerable 
attention by researchers in the computational fluid dynamics area [1-3] and by researchers in the 
experimental fluid dynamics area [4-6]. The problem is of vital importance to the dynamic stability 
and controllability of missiles and fighter aircraft. When flow asymmetry develops, it produces 
side forces, asymmetric lifting forces and corresponding yawing, rolling and pitching moments that 
might be larger than those available by the control system of the vehicle. 

In several recent papers by the present authors [1, 2], the unsteady, thin-layer, compressible 
Navier-Stokes equations have been used to simulate steady and unsteady, asymmetric vortex flows, 
including their passive control, around cones with different cross-sectional shapes. The emphasis of 
these papers was extensive computational studies of the parameters which influence the asymmetric 
flow phenomenon and its passive control. Since the computational cost associated with the solution 
of three-dimensional-flow problems with reasonable flow resolution is very expensive, all the 
computational solutions were obtained using a locally-conical flow assumption. Such an assumption 
reduces the problem solution to that on two conical planes, which are in close proximity of each 
other, and hence it reduces the computational cost by an order of magnitude. Moreover, such 
solutions still provide extensive understanding of the flow physics since one can use very fine grids 
for reasonable flow resolution. In the present paper, we focus on the three-dimensional asymmetric 
flow problem using a very fine grid with high resolution near the solid boundary. 

2. FORMULATION AND COMPUTATIONAL SCHEME HIGHLIGHTS 

The conservative form of the dimensionless, unsteady, compressible, thin-layer Navier-Stokes 
equations in terms of time-independent, body-conformed coordinates are used. The implicit, upwind, 
flux-difference splitting finite-volume scheme is used to solve the unsteady, compressible, thin-layer 
Navier-Stokes equations. The scheme uses the flux-difference splitting scheme of Roe which is 
based on the solution of the approximate Riemann problem. Boundary conditions are explicitly 
implemented. At the plane of geometric symmetry, periodic conditions are used. Frees tream 
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Figure 1. Asymmetric flow solution around a cone of unit length, short-duration side slip. 



Figure 2. Asymmetric flow solution around a cone of unit length, short-duration side slip. 
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Figure 3. Unsteady asymmetric locally-conical flow solution at different time steps within one 
cycle (cylinder axis is a time axis), At = 10 -3 . 


conditions are specified at the inflow boundaries and first-order extrapolation of the flow variables 
is used at the outflow boundaries. The conical shock enclosing the body is captured as pan of the 
solution. On the solid boundary, the no-slip and no-penetration conditions are enforced and the 
normal pressure gradient is set equal to zero. For the temperature, the adiabatic boundary condition 
is enforced at the solid boundary. The initial conditions correspond to the freestream conditions 
with the no-slip and no-penetration conditions on the solid boundary. 

3. COMPUTATIONAL APPLICATIONS AND DISCUSSIONS 

Circular Cone 

A 5° -semi-apex angle circular cone of unit length (cone length is the characteristic length) 
is considered. This is the same circular cone which was considered by the authors in Ref. 1 for 
the locally-conical flow solutions. A three-dimensional grid of 161x81x65 in the wrap around, 
normal and axial directions, respectively, is generated by using a modified Joukowski transformation 
at axial stations. The grid is clustered algebraically in the normal direction of the body using a 
geometric series with minimum grid spacing of 10* at the cone vertex and 10~ 5 at the axial station 
of unit length. The cross-flow grid size of 161x81 is the same grid size which was used for the 
locally-conical flow solutions of Ref. 1. 

With the flow conditions set at a = 20°, = 1.8 and Re = 10 3 , which are the same conditions 

as those of the locally-conical flow of Ref. 1, the three-dimensional solution produces a symmetric 
steady flow, unlike the locally-conical solution which produces asymmetric steady flow. Next, the 
search is directed at obtaining asymmetric flow solutions for the three-dimensional cone flow. In 
Fig. 1, we show the solution in the form of total-pressure loss for the same cone at a = 40°, Moo 
= 1.4 and R e = 4xl0 6 . It is seen that the solution is asymmetric and is nearly self-similar over 
a long axial distance of the cone length. This solution is obtained using a short-duration side-slip 
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Figure 4. Total-pressure-loss contours and surface-pressure coefficient at different axial stations* 
a cone of unit length, a = 40°, = 1.4, R« = 8xl0 6 . 



Figure 5. Asymmetric flow solution around a cone-cylinder configuration 1:1. 
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Figure 6. Surface pressure, total-pressure-loss and Mach contours at different axial stations, 
cone-cylinder configuration 0.5:0.5, a = 46.1°, Moo = 1-6, R« = 6.6 x 10 6 , comparison 
with experimental data (Ref. 5). 


disturbance. When the residual error drops four orders of magnitude, a 2°-side-slip disturbance is 
applied for 100 iteration steps, then it is removed. Thereafter, the pseudo time stepping- is continued 
until the residual error drops again four to five orders of magnitude and a stable asymmetric solution 
is obtained. 

Figure 2 shows the total-pressure-loss solution for the same cone for a higher Reynolds number, 
R« = 8xl0 6 . The asymmetry of the vortex flow becomes much stronger as compared with the 
previous case. The flow asymmetry of this case changes sides along the axial distance and a 
complete wave length of flow asymmetry is formed between the third and ninth cross-flow planes. 
Strong spatially shed vortices exist in the flowfield. This solution is qualitatively similar to the 
unsteady asymmetric locally-conical flow solution at different time steps [1] which is depicted in 
Fig. 3 on a cylinder with the axis of the cylinder representing time. The behavior of the flow 
asymmetry over one period in Fig. 3 is qualitatively similar to the behavior of the flow asymmetry 
QYgj- one wave length in Fig. 2. Figure 4 shows the total-pressure-loss contours and surface pressure 
coefficient at different axial stations for the case of Fig. 2. The solutions at axial stations of X/L = 
0.2 and 0.9 are almost the same (the total pressure losses are drawn to a scale given by the ratio of 
the circular diameters at X/L = 1 station and the local axial station). The flow asymmetry between 
these two stations represents a full wave length. 
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Circular Cone-Cylinder Configurations 

To address the issue of the effect of cylindrical afterbody length on the flow asymmetry a 
cylindrical afterbody of different lengths is added to the unit-length conical forebody. The flow 
around the resulting cone-cylinder configurations is solved with the ^conditions of o - 40 M 
-14 and R, = 4x 10 6 , which are the same flow conditions of the isolated unit-length cone o Fig. 1. 
TOe lengths of the cylindrical afterbody are chosen as 1, 1.5 and 2. The source of flow thsturbance 
ri the same shori-dunttion 2” -side-slip disturbance. For the cooe-cyltnder coitiiguranon of 1:1 (cone 
length' cylinder lengtit). Fig. 5 shows a very strong asymmetric flow on the cone^n comparison 
with the flow asymmetry of the isolated cone of Fig. 1, and on the cylindrical afterbody as well. It 
Tould be noted that inride the conical shock smrounding the cone-cyhnder configuration. subsonrc 
"^and hence the downstieaa cyltndncal-afterbody boundary has anupstrcam effect 

The cylindrical afterbody has dual effecB which in^^» «^m«n-. ^ 

Tc sssrflrjsr =-/-? si = -csr/is 

asymmetry. For the cone-cylinder configurations of 1.1.5 and 1.2, th ymm try 

the flow becomes unsteady [Ref. 3]. data For this 

Next we show a comparison of the computed results with available expenme 

prnpo e we consider the cone-cyiinder configuration of 0.5. 0.5 which was «P^»^ 

P ium>. The configuration angle of atiack is 46.1«, dm Mach numbm n U 

number based on the total configuranon length <c0 * 161x81x65 Figure 6 shows the 

anex anele is 9.5°. The problem is solved using a gnd size of 161 x six . g 

surface-pressute coefficient altmg with the 
leisure coefficient arc in good agreement on all the axial 

rides in the downstream direction as it is shown by the results of aroal sunons at 0.475 and 0.775. 
This comparison conclusively validates our computed results and the gnd size. 

4 - eve rri* hnpomm t^issue^ arc ^addre ssed in tire present sn^y. By incre^ng. fieflow Refolds 

SS- : i 

« By ling a cylindrical afterbody re tire conical forcbody, tire flow asymm.tr, 

grid used are conclusively validated. 

fe y NASA Research Centcr under 

1 _ 994 . The computations have been carried out using the CRAY 2 (Voyger) of the N g ey 

Research Center. 
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Prediction of Steady and Unsteady Asymmetric 
Vortical Flows Around Circular Cones 
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Steady and unsteady, supersonic asymmetric vortical Hows and their passive control around circular cones are 
considered in this paper. These problems are formulated by using the unsteady, compressible, single and double, 
thtn-layer Navier- s tokes equations. The equations are solved by using an implicit, upwind, flux-difference 
sp ting, inile- volume scheme, either in a pseudotime stepping or in an accurate time stepping. An implicit, 
approximately factored, central-difference, finite-volume scheme has also been used to validate some applica- 
tions of the upwind scheme. Local conical flows are assumed for the computational applications presented in 
is paper. Steady asymmetric vortical flows have been predicted by using random and controlled disturbances, 
nsteady asymmetric vortex-shedding flows have also been predicted, for the first time, using time-accurate 
solutions with two different computational schemes. Control of flow asymmetry has been demonstrated 
computationally by inserting a vertical fin in the leeward plane of geometric symmetry. 


Introduction 

I N the high angle of attack (AOA) range, the separated 
vortical flow from forebodies of missiles and fighter air- 
craft may become asymmetric, producing large abrupt 
changes in force and moment coefficients. These abrupt 
changes may exceed the available controllability and lead to 
missile and aircraft spin. Experimental studies of several re- 
searchers 11 have identified four distinct flow patterns about 
slender bodies through a wide AOA range and zero-degree 
side slip. The first pattern develops in the very small AOA 
range, where the flow is attached and the axial flow is domi- 
nant. In the intermediate AOA range, the crossflow becomes 
of the same order of magnitude as that of the axial flow, the 
flow separates on the leeward side, and a symmetric vortex 
pair is formed. As the AOA reaches a high range, the symmet- 
ric vortex pair becomes asymmetric, and the flows stay steady. 
For this asymmetric vortex-flow pattern to occur, it is not a 
necessary condition to have asymmetric separation lines on the 
leeward side of the body. The fourth flow pattern develops at 
a very high AOA range, where asymmetric time-dependent 
vortex shedding occurs either randomly or periodically, simi- 
lar to the von Karman vortex street in two-dimensional* flows 
around cylinders. 

For isolated pointed forebodies, the onset of vortical flow 
asymmetry occurs when the relative incidence (ratio of AOA 
to nose semi-apex angle) exceeds a certain value; e.g., for a 
pointed circular cone, the relative incidence must be higher 
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than two. However, the relative incidence value is not the only 
determinable parameter for the onset of vortical flow asym- 
metry. The onset of vortical flow asymmetry is also a function 
of the freestream Mach number and Reynolds number and the 
shape of the body cross-sectional area as well. Asymmetric 
vortical flow and vortex shedding have also been documented 
for delta wings 1213 at very high relative incidences and low 
subsonic regimes. 

For the critical values of the relative incidence, Mach num- 
ber and Reynolds number, and the shape of cross-sectional 
area, the symmetric flow is unstable. Any small flow distur- 
bance in the form of a transient side slip, acoustic disturbance, 
or similar source of disturbance causes flow instability that 
produces, depending on the flow conditions, either a steady 
asymmetric vortical flow or an unsteady asymmetric flow with 
vortex shedding. In this paper, we present an extensive compu- 
tational study of the steady asymmetric vortical flow and 
unsteady asymmetric flow with vortex shedding to address 
some of the influential parameters as the relative incidence 
and Mach number. 

As the experimental work shows, the mechanisms that lead 
to asymmetric vortex wake are not well understood. However, 
two mechanisms have been established for explaining the evo- 
lution of now asymmetry. 5 - 6 - 9 - 10 The first mechanism applies 
to both laminar and fully turbulent flows. It suggests that flow 
asymmetry occurs due to instability of the velocity profiles in 
the vicinity of the enclosing saddle point that exists in the 
crossflow planes above the body primary vortices. 2 ’ 710 The 
second mechanism suggests that flow asymmetry occurs due to 
asymmetric transition of the boundary-layer flow either at the 
nose in the axial direction or on both sides of the body in the 
crossflow planes. For pointed slender bodies, the first mecha- 
nism produces higher side forces than those produced by the 
second mechanism. These results have conclusively been 
shown through the experimental work of Lamont 8 9 on 2-diam 
and 3.5-diam tangent ogive noses with cylindrical afterbody. 
An extensive review of the steady and unsteady vortex-induced 
asymmetric loads is given by Ericsson and Reding in Ref. 11. 

Several attempts have been carried out to computationally 
simulate asymmetric vortical flows around slender bodies of 
revolution. Early computational work on conical flows has 
been published in Refs. 14 and 15. Graham and Hankey 16 
presented the first three-dimensional Navier-Stokes computa- 
tions for asymmetric flow around a cone-cylinder bodv ar 30- 
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deg angle of attack, 1.6 freestream Mach number, and 
0.4 x 10 6 Reynolds number. The MacCormack explicit finite 
difference scheme was used for the computations on a rela- 
tively coarse grid of 26 x 30 x 60. A very small perturbation is 
induced by the truncation error of finite difference algorithm 
that triggers an instability of the saddle point above the body 
(first mechanism for asymmetry). Hence, the instability is 
induced by numerical bias that is physically amplified to pro- 
duce flow asymmetry. By switching the order of spatial differ- 
encing in the predictor and corrector sweeps, the asymmetry 
was reversed. 

Degani and Schiff 17 used the thin-layer, Reynolds-averaged, 
Navier-Stokes equations to compute asymmetric vortical flow 
around an ogive-cylinder body. They found that flow asym- 
metry can be obtained by introducing an asymmetric distur- 
bance very close to the body nose. The disturbance they used 
was in the form of a small jet that was blown from one side of 
the body near the nose. However, when the jet was turned off, 
the numerical solution unfortunately showed that the flow 
recovered its symmetry. The authors of the present paper 
believe that the problem is attributed to the smallest scale of 
the grid at the solid boundary and the damping effect of the 
numerical dissipation in the axial direction, in addition to the 
grid-fineness distribution. 

Marconi 18 used the Euler equations to solve for supersonic 
flow past a circular cone in conjunction with a “forced separa- 
tion model,” which was used by Dyer, et al. 19 The pseudotime 
stepping was carried out until the residual error reached ma- 
chine zero while the flow was symmetric. Proceeding with the 
time stepping, vortex-flow asymmetry was obtained and 
stayed stable thereafter. It is believed that the asymmetry was 
triggered by the machine round-off error, which acted as a 
disturbance to the saddle point in the flowfield. In a later 
paper, Siclari and Marconi 20 used the full Navier-Stokes equa- 
tions to solve for supersonic asymmetric flows around a 5-deg 
semiapex angle cone over a wide range of angles of attack. 

Very recently, Stahl 21 conducted experimental studies of the 
low-speed flow around a circular cone of 8-deg semiapex angle 
circular cone in the angle of attack range of 15-50 deg at a 
Reynolds number of 7800 based on the base diameter. The 
onset of flow asymmetry was observed at 35-deg angle of 
attack. He has shown that the flow asymmetry can be sup- 
pressed by inserting a fin along the leeward plane of geometric 
symmetry with its edge along a ray through the apex. The 
minimum fin height for this purpose was found to be equal to 
the local radius of the cone. 

In this paper, the supersonic, steady and unsteady, asym- 
metric vortical flows around circular cones are studied using 
the unsteady, compressible, single thin-layer, Navier-Stokes 
equations. Two computational schemes are used to solve the 
equations. The first, which is the main scheme used in this 
paper, is an implicit, upwind, flux-difference splitting, finite- 
volume scheme. The second, which is used to validate certain 
cases of the upwind scheme, is an implicit, approximately 
factored, central-difference, finite-volume scheme. Pseudo- 
time stepping is used for steady flows and time-accurate step- 
ping is used for unsteady flows. Some of the influential 
parameters for flow asymmetry, such as the relative incidence 
and Mach number, are addressed. A flow case of passive 
control of flow asymmetry is also studied using the unsteady, 
compressible, double thin-layer, Navier-Stokes equations. 


where 


Q = -j = -j [p,pu u pu2,pui,pe\ (2) 

E m m inviscid flux 


= (\/J)[pU m ,pU x U m + 

+ {f>e + p)U n ,Y, m = 1,2,3 (3) 

(E v ) 2 = viscous and head-conduction flux in direction 

= (l//)[0,d*£ 2 7*|, d*£ 2 T*2, d*£ 2 7 * 31 (w„7 kn~Qk)\ W 

(£ v )j = viscous and head-conduction flux in £ 3 direction 

= (l/y)[0,d* { 3 7*1, d*£ 3 7*2» d*£ 3 7*3, d*£ 5 {Un T kn~ Qk)] (^) 

u m = d k ruk * 6 > 

The first element of the three momentum elements of Eq. (5) 
is given by 

where 

<t> = a k ea k e. <*> 


The second and third elements of the momentum elements are 
obtained by replacing the subscript 1, everywhere in Eq. (7), 
with 2 and 3, respectively. The last element of Eq. (5) is given 
by 

a k e (u„T kn -q k ) - ^ + *[ ‘/’apW + " 2 + “3) 


+ 


1 

(7 - 1 )Pr 


d(g 2 ) ] 

Be JJ 


(9) 


where 


w=d n eu n < 10 > 

For Eq. (4), in the case of double thin-layer, Navier-Stokes 
equations, the elements are given by equations similar to Eqs. 
(7-10) with the exception of replacing £ 3 by £ . The double 
thin-layer, Navier-Stokes equations are used only for the pas- 
sive control of flow asymmetry since the existence of the fin 
creates a second thin layer that is perpendicular to the cone 
thin layer. The reference parameters for the dimensionless 
form of the equations are £, L/a *, p«, and p* for the 
length, velocity, time, density, and molecular viscosity, re- 
spectively. The Reynolds number is defined as Re = f>„V„L/ 
n mt and the pressure p is related to the total energy per unit 
mass and density by the gas equation 

p = (-y — l)p [e— yi(U\ + u* + W3 2 )] (H) 


Formulation 

The three-dimensional compressible viscous flow around 
the body is governed by the conservative form of the dimen- 
sionless, unsteady, compressible, double thin-layer, Navier- 
Stokes equations. In terms of time-independent, body-con- 
formed coordinates £', J 2 , and £ 3 , the equations are given by 

dQ 3E> B(E v h Bfoh _ 0 5= 1,2,3, (1) 

di + dv se Be 


The viscosity is calculated from the Sutherland law 

C= 0.4317 (12) 

and the Prandtl number P, = 0.72. 

In Eqs. (1-10), the indicial notation is used for convenience. 
Hence, the subscript k and n are summation indices, the 
superscript or subscript s is a summation index, and the super- 
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script or subscript m is a free index. The range of k , n , s, and 

m is 1-3, and d k * (d/d Xt ). ..... . . T . in 

Boundary conditions are explicitly implemented. They in- 
clude inflow-outflow conditions and solid-boundary condi- 
tions. At the plane of geometric symmetry, periodic condi- 
tions are used for symmetric or asymmetric How applications 
on the whole computational domain (right and left domains). 
At the far-field inflow boundaries, freestream conditions are 
specified since we are dealing with supersonic flows, whereas 
at the far-field outflow boundaries, first-order extrapolation 
from the interior points is used. On the solid boundary, the 
no-slip and no-penetration conditions are enforced; u,=u 2 - 
Ui = 0, and the normal pressure gradient is set equal to zero. 
For the temperature, the adiabatic boundary condition is en- 
forced on the solid boundary. The initial conditions corre- 
spond to the uniform flow with u, = u 2 = = 0 on the solid 

boundary. . _ • 

For the passive control applications using a vertical tin in 
the leeward plane of geometric symmetry, solid-boundary 
conditions are enforced on both sides of the fin. 


Highlights of Computational Schemes 

The first computational scheme used to solve the unsteady 
compressible, single or double thin-layer, Navier-Stokes equa- 


tions is based on the Roe inviscid flux-difference splitting 
scheme. In this scheme, the Jacobian matrices of the inviscid 
fluxes, A s =(dE/dq), s= 1-3, are split into left and right 
fluxes according to the signs of the eigenvalues of the inviscid 
Jacobian matrices. Flux limiters are used to dampen the nu- 
merical oscillations in regions of large changes of the gradients 
of the flowfield vector. The viscous and heat transfer terms 
are centrally differenced. The resulting equation is solved by 
using approximate factorization in the £ l , and £ direc- 
tions. The computational scheme is coded in the computer 

program CFL3D. ... * 

The second computational scheme is an implicit, approxi- 
mately factored, centrally differenced, finite-volume scheme. 
Added second-order and fourth-order dissipation terms are 
used in the difference equation on its right-hand side terms, 
which represent the explicit part of the scheme. The Jacobian 
matrices of the implicit operator on the left-hand side of the 
difference equation are centrally differenced x n space, an 
implicit second-order dissipation terms are added for t e 
scheme stability. The left-hand side operator is approximately 
factored, and the difference equation is solved in three sweeps 
in the £', £ 2 , and £ 3 directions, respectively. The computa- 
tional scheme is coded in the computer program 1CF3U. l ne 
ICF3D code is used to verify some of the applications of the 



— CFL3D 
— ICF3D 




Fig. i 


Symmetric flow 


solutions for a circular cone. 


o = 10 deg, M. = 1.8. Re « 10 5 (validation case). 


Fig. 2 



Symmetric 


— CFL3D 





Steady asymmetric flow solutions for a circular cone due to random 


disturbances. «r = 20 deg, = 1.8. R, = 10 s (validation case). 
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CFL3D code; namely the cases of Figs. 1 and 2. For the 
problem of passive control of now asymmetry, the double 
thin-layer, Navier-Stokes equations have been solved using the 

CFL3D code. , . 

Since the applications in this paper cover local-conical 1 lows 
only, the three-dimensional scheme is used to solve for locally 
conical nows. This is achieved by forcing the conserved com- 
ponents of the How vector field to be equal at two planes of 
x = 0.95 and 1.0. The validity of local-conical-now assump- 
tion is discussed in the next section. 

Validity of the Local-Conical-Flow Assumption 

The solutions presented in this paper are called local-conical 
solutions, which are obtained by equating the conserved com- 
ponents of the nowfield vector, in the three-dimensional 
scheme, on two crossnow planes that are in close proximity to 
each other at a selected location. Once this location is specified 
(x = 1 0 in the present applications), the now Reynolds num- 
ber is determined and the time scale, for time-accurate solu- 
tions, is also determined. The resulting solution is a local-con- 
ical solution at the specified location. It is not a global-conical 
solution. The locally conical equations can be shown by con- 
sidering the conservative form of the Navier-Stokes equations 
in the Cartesian system 

dq + djE-E J, = Q' ,• = 1-3 (13) 

dt d x, 

By introducing the conical coordinates 


coordinates). Hence, Eq. (18) is not self-similar, and there- 
fore it does not represent a global-conical flow. However, if 17* 
is set equal to a constant c, then one can consider Eq. (18) to 
represent a local-conical flow around 771 = c. The resulting 
solution using Eq. (18) with V} = c represents a local-conical 
solution with a Reynolds number and a time that are scaled by 
the constant c. It should be noted that if the flow is steady and 
inviscid, then Eq. (18) becomes self-similar, and hence it rep- 
resents a global-conical How. In the present paper, we indi- 
rectly solve Eq. (18) at a fixed location of unity. This is 
achieved in the three-dimensional now equation, Eq. (1), by 
equating the elements of the nowfield vector at two planes in 
close proximity to each other. In this paper, w-e selected these 
planes to be located at * - 0.95 and 1.0. In other numerical 
experiments, we use the plane locations at x = 0.995 and 1.0. 
The results of these experiments were in excellent agreement 
with those of the present paper. 

Computational Studies 

Supersonic fiows about a 5 deg semiapex angle circular cone 
at a Reynolds number of 10 5 have been considered. A grid of 
161 x 81 points in the circumferential and normal directions is 
used throughout the present applications. The grid is gener- 
ated by using a modified Joukowski transformation with a 
geometric series for the grid clustering near the solid 
boundary. The minimum grid length is 10“ 4 at the solid 
boundary, and the maximum radius of the computational 
domain is 21 r, where r is the radius of the circular cone at the 
axial station of unity. 


V 1 


£1 



nl = XiXi 


(14) 


and using the chain rule to express Eq. (13) in terms of the 
conical coordinates, we get 


*£ _ Erh + -L(£-£ v)2 

6 dt dr) t <7772 

+ ^ (E-E v h + 2(7— A) = 0 

0 - dr) } 


(15) 


where 


6 - v 1 +777 + 172 

£,=£,- T? 1 £3, £: = £:- V2E3 

£3 = £3 + t ? i£i + V 2^2 


?«£, < 16 > 

E v \ = £ V 1 - T?l£v3 

£ v2 = E\2 — l?2^v3 

E vi = £V 3 + T?1 £1-1 + 

/, = £„ 3 < 17 > 

The conical flow condition requires that the flow variables be 
independent of the coordinate 773. If this condition is imposed 
in Eq. (15), by dropping the derivatives with respect to tj,, the 
equation reduces to 

+ £(£-£,,), + -^-(E-Evh + 2(7 -E) = 0 (18) 

6 St dr) ] dr) 2 

It is clearly seen that Eq. (18) still has m dependence in the 
unsteady term and the viscous and heat-flux terms (one can see 
the explicit dependence of the viscous and heat-flux terms on 
773 by transforming the elements of these vectors to the conical 


Steady Symmetric Flows 

Figure 1 shows steady symmetric vortical-flow solutions for 
the circular cone at 10 deg angle of attack and 1.8 freestream 
Mach number. In the figure, we show comparisons of the 
results of the CFL3D and 1CF3D codes. The results include 
the residual error versus the number of iterations, the cross- 
flow velocity, the total-pressure-loss contours, and the sur- 
face-pressure coefficients. It should be noted here that the 
angle 6 in the C p figure is measured from the leeward plane of 
geometric symmetry in the clockwise direction. The agreement 
of the results of the two code is excellent, and the results are 
in full agreement with those of Siclari and Marconi. 1 

Steady Asymmetric Flow 

Round-Off and Truncation Error Disturbances 
The cone angle of attack is increased to 20 deg while all the 
other flow conditions are kept fixed. Figure 2 shows the results 
of the CFL3D and 1CF3D codes. In the residual error figure, 
the CFL3D code shows that the residual error drops 10 orders 
of magnitude within 2500 iteration steps. Thereafter, the error 
increases by six orders of magnitude. The flow is symmetric 
during this 5000 iteration steps. Next, the error drops down by 
another six orders of magnitude and stays constant for 2500 
iteration steps. The now becomes asymmetric and stable. The 
1CF3D code shows that the residual error drops five orders of 
magnitude in the first 3000 iteration steps, increases two or- 
ders of magnitude in the next 2000 iteration steps, and then 
drops down by three orders of magnitude within the next 5000 
iterations. The flow solution goes through a symmetric un- 
stable solution and then to an asymmetric stable solution. The 
pressure-coefficient figure for the two codes is the same over 
the full range of the circumferential angle 8. The suction 
pressure in the range of 6 = 0-90 deg is lower than that of the 
range of 8 = 270-360 deg. The crossflow velocity and total- 
pressure-loss contours for the two codes are also in excellent 
agreement. They show the nature of the now asymmetry and 
its details. The results are in complete agreement with those of 

Ref. 20. J . „ 

Since the residual error of the CFL3D code is much smaller 
than that of the ICF3D code after the first 2500 iterations, the 
disturbance that triggered the asymmetry in the first code is 
attributed to the machine round-off error, while the distur- 
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Fig. 3 



Steady asymmetric flow solutions for a circular cone due to 


±2 deg transient side slip, a = 20 deg, - 1-8* Re - 


10 s . 


Fig. 4 



Steady asymmetric flow solutions for a circular cone due to ± 5 deg transient side slip, a = 20 deg. - 1.8. Re - 10 s 


bance that triggered the asymmetry in the second code is 
attributed to the truncation error of the scheme (since there is 
a bias due to the spatial marching direction). Both distur- 
bances are random in nature. However, irrespective of the 
source of disturbance, the final asymmetric stable solution is 
the same. 

Coni rolled Transient Side-Slip Disturbances 

In Figs. 3 and 4, we show steady asymmetric flow solutions 
due to transient side-slip disturbances of ± 2 and ± 0.5 deg. 
The residual-error figures show a drop of seven orders of 
magnitude in the first 2000 iterations. At this step, a side-slip 
disturbance is imposed for six iteration steps, then it is re- 
moved. Irrespective of the magnitude or the sign of the side- 
slip disturbance, the residual error increases by six orders of 
magnitude, then it drops down very rapidly. A stable asym- 
metric flow solution is obtained. The asymmetric solutions 
corresponding to the ± 2 deg side-slip disturbances are mirror 
images of each other, as can be seen from the figures of the 
surface-pressure coefficient, crossflow velocity, and total- 
pressure-loss contours. The corresponding asymmetric solu- 
tions with the ±0.5 deg side-slip disturbances are exactly the 
same as those of the ± 2 deg side-slip disturbances. Moreover , 


the Final asymmetric solutions of the ± 2 deg and ± 0.5 deg 
side-slip disturbances are the same as those of Fig. 2. 

Again, this numerical experiment shows that the same phys- 
ical flow asymmetry is obtained. 

Unsteady Asymmetric Vortex Shedding 

In the present case, the angle of attack is increased to 30 deg 
and all the other flow conditions are kept the same as those of 
the cases above. Figure 5 shows the results of this case. 

Here, we show the history of the residual error and the lift 
coefficient up to the 15,700 time step. First, pseudo-time 
stepping was used up to 10,000 iterations, and the solution was 
monitored every 500 iterations. The solution showed that the 
asymmetry was changing from the left side to the right side, 
which indicated a possibility of unsteady asymmetric vortex 
shedding. The residual error was also oscillating. The compu- 
tations were repeated starting from the 3,500 iteration step 
using time-accurate calculations with At — 10 3 . The residual- 
error and lift-coefficient figures show the time history of the 
solution. It is seen that the residual error and the lift coeffi- 
cient show a transient response that is followed by a periodic 
response. Figure 5 shows also snapshots of the time history of 
the solution for the total-pressure-loss contours and surface- 
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15,700 



Fig. 5 Unsteady asymmetric How solution (using flux-difference splitting) with vortex shedding for a circular cone during periodic flow responses; 
a = 30 deg, M*, = 1.8, R, = 10 s , A/ = 10' 3 . 


pressure coefficient. The solutions are shown every 100 time 
steps starting from the time step of 15,000. At n = 15,000, the 
asymmetric flow is seen with an already shed vortex from the 
right side. As time passes, the shed vortex is convected in the 
flow and the primary vortex on the left side stretches upwards 
while the primary vortex on the right gets stronger, as it is seen 
from the surface pressure figures. At n = 15,600, the primary 
vortex on the left side is about to be shed. At n = 15,700, the 
primary vortex on the left side is shed in the flowfield. It 
should be noticed that the solution at n - 15,700 is exactly a 
mirror image to that at n = 15,000. The solution from 15,000- 
15,700 represents the first one-half the cycle of shedding. The 
solution from 15,700-16,400 (not shown) represents the sec- 
ond one-half the cycle. The periodicity of the shedding motion 
is conclusively captured. The period of oscillations is 10 ~ 3 x 
1,400 steps = 1.4 that produces a shedding frequency of 4.400 
(Strouhal number). This solution is obtained by using the 
flux-difference splitting (FDS) scheme. 

Very recently, a researcher in the computational simulation 
area of asymmetric flows claimed that he had applied the 
flux-vector splitting (FVS) scheme of the CFL3D code to the 
present flow case. His solution showed that the flow was 


steady and symmetric. A statement of his results was commu- 
nicated to us and we were asked to respond. Therefore, we 
recomputed the present flow case using the FVS scheme of the 
same CFL3D code. In Fig. 6, we show the results of the 
time-accurate solutions using the FVS scheme using the same 
grid. Using the FVS scheme, the flux limiters were turned on, 
and as can be seen from the logarithmic-residual curve, the 
solution becomes symmetric and steady after 5000 time steps. 
Next, the flux limiters are turned off, and the solution shows 
a transient response up to 12,000 time steps. Thereafter, the 
solution becomes periodic with periodic asymmetric vortex 
shedding. The solution was monitored every 100 time steps, 
and the results from n = 13,900-14,600 are shown. Although 
the process of adjusting the time instants is difficult to match 
those of the FDS solution, it is seen that the captured snap- 
shots of the FVS solution almost match those of the FDS 
solution. Comparing the FVS solutions at n = 13,900 and 
14,600, it is seen that they are mirror images of each other. 
Hence, periodic flow response has been achieved with a period 
of 1400 x 10 “ 3 = 1.4, which is exactly the same period of 
shedding as that of the FDS solution. This pinpoints the high 
numerical dissipation effect of the FVS scheme when the flux 
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limiters are turned on. The resulting numerical dissipation in 
the FVS is large enough to dampen the random disturbances 
of the flow solution. By turning off the flux limiters in the 
FVS scheme, the random disturbances can grow, producing 
the asymmetric unsteady vortex shedding. This also shows 
that the FDS scheme, even with the flux limiters turned on, is 
less dissipative than the FVS scheme. These results conclu- 
sively explain the erroneous claim of steady flow made by the 
previously mentioned researcher. 

Steady Asymmetric Flow at Different Mach Numbers (Effect of M„) 
Figure 7 shows the effect of the freestream Mach number 
(M„ = 2.2, 2.6, and 3.0) on the convergence history, surface 
pressure! crossnow velocity, and total-pressure-loss contours 
for the circular cone at 20 deg angle of attack. At A/» = 2.2, 
the residual error shows that the stable asymmetric flow is 
obtained within the same number of iterations as that of the 
= 1.8 case. At = 2.6, the residual error shows that the 
stable asymmetric flow is obtained after a large number of 
iterations. And at = 3.0, no asymmetric flow was cap- 
tured, the flow stayed symmetrically stable. The surface pres- 




iterations 


sure figures show that the asymmetry gets weaker as the Mach 
number is increased. This conclusion is clearly seen from the 
crossflow velocity and the total-pressure-loss figures. It should 
be noted that since the nature of disturbance is random, flow 
asymmetry changes sides as the Mach number increases until it 
disappears. 


Passive Control of Flow Asymmetry 

Figure 8 shows the passive control of flow asymmetry by 
inserting a vertical fin in the leeward plane of geometric sym- 
metry. The fin height is equal to the cone local radius r . Here, 
the double thin-layer, Navier-Stokes equations are used to 
obtain these results. The flow Mach number is kept at 1 .8 and 
the angle of attack is 20 deg. The How is completely symmetric 
as can be seen from the figures of the surface-pressure coeffi- 
cient, total-pressure-loss contours, and crossflow velocity. A 
blow-up of the cross-flow velocity at the fin-cone juncture 
shows two corner recirculating bubbles of exactly the same 
size. This case has been obtained after 24,000 iteration steps. 
Again, this is the first time such a computational simulation of 




Fig. 6 Unsteady asymmetric flow solution (using finx-vector splitting) with vortex shedding for a circular cone during periodic flow response: 
a = 30 deg, Mao — 1 .8, Re = 10 s , A/ = 10 
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Fig. * Control of .symmetric flow of . circular cone using a vertical fin in the leew.nl geometric plane of symmetry, a = 20 deg. - 1.8 
R, - 10 s , h = r <r ■ radius of circular section). 


the passive control of the flow asymmetry has been presented. 
The results are in full agreement with Stahl s experimental 

study. 21 

Concluding Remarks 

This paper presents extensive computational study and sim- 
ulation of steady and unsteady asymmetric vortex flow around 
circular cones. A systematic study has been carried out to 
show the effects of angle of attack and Mach number. The 
study shows that the flow asymmetry is independent of the 
type or level of the disturbance. For the controlled transient 
side-slip disturbance, the solution is unique. For the uncon- 
trolled random disturbance, the solution is also unique with 
the exception of having the same asymmetry changing sides on 
the cone. It conclusively shows that periodic vortex shedding 
has been captured at larger angles of attack. The unsteady 
asymmetric vortex-shedding solution has been substantiated 
by using two different computational schemes. It also shows 
that as the Mach number increases, the vortex flow asymmetry 
gets weaker until it disappears. The possibility of passive 
control of flow asymmetry has also been demonstrated. Many 
of the cases presented here are obtained for the first time, in 
particular, the asymmetric vortex shedding cases and the cases 
of passive control of flow asymmetry. 
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